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■ Abstract 

_ We introduce the GENJI program (Gamma-ray Emitting Notable AGN Monitoring by Japanese VLBI), which is 

f~i . a monitoring program of gamma-ray bright AGNs with the VERA array (VLBI Exploration of Radio Astrometry). 

Pj The GENJI programme aims a dense monitoring at 22 GHz towards the 7-ray emitting active galactic nuclei (AGNs) 

Q . to investigate the radio time variation of the core and possible ejection of new radio component, motion of jets, and 

^ ■ their relation with the emission at other wavelengths especially in 7-rays. Currently we are monitoring 8 7-ray- 

! emitting notable AGNs (DA 55, 3C 84, M 87, PKS 1510-089, DA 406, NRAO 530, BL Lac, 3C 454.3) about once 

■ every two weeks. This programme is promising to trace the trend of radio time variation on shorter timescale than 
I conventional VLBI monitoring programme and to provide complimentary data with them (e.g., MOJAVE, Boston 

■ University Blazar Project). In particular, we successfully coordinated quick follow-up observations after the GeV 
^ ' 7-ray flare in NRAO 530 and 3C 454.3 reported by the Fermi Gamma-ray Space Telescope. Here we present the 

. initial results of morphology and light curves for the first 7-month operation. 

• Key words: galaxies: active, galaxies: jets, galaxies: individual (DA 55, 3C 84/NGC 1275, M 87, PKS 1510- 

089, DA 406, NRAO 530, BL Lac, 3C 454.3), radio continuum: galaxies 



o: 

Introduction At radio wavelengths, the high angular resolution provided 

■ by Very Long Baseline Interferometry (VLBI) is essential to 
The non-thermal emission from the relativistic jets emanat- identify the radio counterpart of the 7-ray emitting site in AGN 

ing from the super massive black holes usually dominates at a jets. In conventional framework of spectral energy distribu- 

• '~j ■ wide range of electromagnetic spectrum from radio to 7-rays. tion (SED) modeling with one-zone synchrotron-self Compton 

rN '. Some emissions appear up to very high energy range even at and/or external Compton model, the blazars should have a 

■ GeV and TeV 7-ray energy ranges. The location of high energy subpc-sized emission region or even smaller size (Kubo et al. 
- ■ ' emission and its production mechanism are long standing prob- 1998) and the emission region is possibly associated with the 

lem in AGN jet physics. With the recent progress by the Large region close to the jet base. Ultimately, VLBI has a potential 

Area Telescope on-board Fermi Gamma-ray Space Telescope to identify the exact location of 7-ray emitting region. In ad- 

and new generation Cherenkov telescopes such as H.E.S.S., dition, VLBI is a unique tool to measure the jet motion, which 

MAGIC, VERITAS, we have a new opportunity to explore the allows us to investigate whether the 7-rays beamed with the 

jet physics in connection with 7-ray emission. The number same Lorentz factor as indicated by the VLBI motion, 

of AGN-hosted 7-ray sources have drastically increased com- The 7-ray luminosity shows a good correlation with the ra- 

pared to the era of the EGRET/CGRO (Nolan et al. 2012), dio luminosity over more than 4 orders of magnitude (Lister et 

including new discoveries of 7-ray emission from misaligned al. 2011). This seems to indicate a tight connection between 

AGNs (e.g., Abdo et al. 2010a). While a number of compet- 7-ray and radio properties in AGN jets. So far considerable 

ing models have been proposed (e.g., Sikora et al. 2009 and efforts have been expended to do intensive VLBI monitoring 

references therein), the 7-ray emission mechanism and the lo- in the context of radio-7-ray connection. However, not all the 

cation of emitting region are still disputed. Quasi-simultaneous 7-ray events in blazars and radio galaxies have a clear radio 

multiwavelength study is a key to discriminate such models. counterpart as in the case of the GeV- TeV 7-ray flare in 3C 

279 (Abdo et al. 2010b) and GeV 7-ray flare in 3C 84 (Nagai 
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some 7-ray flares seem to be related to changes in radio band 
(e.g., Marscher et al. 2010a; Orienti et al. 2011) while others 
show no relation (e.g., D'Ammando et al. 2009). Therefore, 
the location of radio counterpart of the 7-ray emission is still 
an open question and additional validations in more samples 
are highly awaited. 

Gamma-ray Emitting Notable AGN Monitoring by Japanese 
VLBI (GENJI) programme aims dense monitoring of the 7-ray 
bright AGNs using the VLBI Exploration Radio Astrometry 
(VERA). VERA aims at Galactic maser astrometry, but for 
technical reason VERA needs to observe a bright calibrator 
every ^ 80 min. We use this calibrator time as the slot for 
GENJI sources. Since most of the VERA projects are car- 
ried out at 22 GHz, GENJI observation is also basically lim- 
ited to 22-GHz observation. But, sometimes we can perform 
43-GHz observation depending on the aim of VERA project 
(e.g., SiO maser observation) ' . The goal of GENJI program is 
to study the correlation of time variability between 7-ray and 
radio components including the radio core and some other jet 
components in order to identify a possible radio counterpart 
of the 7-ray emission. We also aim to measure the apparent 
motion of the jet allowing us to calculate the jet speed. We 
focus on the time variation on the timescale of shorter than 1 
month and quick follow-up after the 7-ray flare. Currently, we 
are carrying out this monitoring once two-weeks. Such a dense 
monitoring is a unique at 22 GHz in the northern hemisphere 
{cf. TANAMI in the southern hemisphere: Ojha et al. 2010) so 
that we can obtain complementary data with other monitoring 
projects (e.g., MOJAVE^: Lister et al. 2009; Boston University 
Blazar Projects'*). 

The GENJI observations started in 2010 November. We have 
already detected radio flaring associated with a GeV flare in 
PKS 1510-089 (Orienti et al. 2012). In addition, we have also 
successfully coordinated observations right before and after the 
GeV flare in 3C 454.3 and NRAO 530, and these results are 
now in preparation. In this paper, we present a quality of the 
data and initial results for first 7-month data. Throughout the 
paper, we adopt the following cosmological parameters; = 
70.2 km sec^i Mpc^^, VIm = 0.27, and f^A = 0.73 (Komatsu 
et al. 2011). 

2. Sample 

We do not demand statistical completeness of AGN jets, 
and we rather aim dense sampling for particular sources. The 
GENJI sample consist of seven bright blazars and two radio 
galaxies (3C 84 and M 87). The sample is selected from the 
VLBA 2-cm Survey (Kellermann et al. 1998) and MOJAVE 
programme (Lister et al. 2009), whose expected correlated flux 
density is brighter than 1 Jy to fulfill a role of fringe finder or 
delay calibrator for the VERA project observations, identified 
as GeV 7-ray source by Fenni/LAI. Because GENJI observa- 
tions should be carried out within a short time slots not to dis- 
rupt total observing time for the VERA project by antenna slew 
as much as possible, we can basically choose the sources close 

' The VERA has no frequency agihty, and therefore we cannot do 22 and 

43 GHz observations at the same time. 
^ http: / /www. physics. purdue.edu /MOJAVE / index.html 
' http://www.bu.edu/blazars/research.html 



to the target source of the VERA project on the celestial plane. 
We finally selected 8 source listed in Table 1 . The plots of typi- 
cal interferometric uw-coverage are shown in Figure 1 . OJ 287 
has been included in the GENJI sample since 201 1 December. 

3. Observation and Data Reduction 

Observations were carried out using the VERA 4 stations. 
The monitoring was suspended from 2011 March 11 to 2011 
April 19 because of a big earthquake hitting north-east Japan. 
From 2011 June to 2011 August, the VERA was not able to 
be used because of the maintenance. Most of the sources were 
observed at 22 GHz. Sometimes 43-GHz observations were 
also performed, but we do not include the 43-GHz result in 
this paper On-source time for each observation was typically 
30 minutes, spreading over 4-6 scans at different hour angle. 
Each source of the sample was observed almost once every two 
weeks. 

Left hand circular polarization (LHCP) was received and 
sampled with 2-bit quantization, and filtered using the VERA 
digital filter unit (Iguchi et al. 2005). The data were recorded at 
a rate of 1024 Mbps. The baseband allocation depends on the 
observation, but mostly providing a band-width of 256 MHz in 
which 14 IF-channels per a total of 15 IF-channels of 16 MHz 
bandwidth were assigned to the GENJI source. Because of a 
significant signal loss (^ 30%) due to the analog filter imple- 
mented in front of the digital filter, we do not use the data of 
12th, 13th and 14th IFs. Correlation processes were performed 
with the Mitaka FX correlator (Chikada et al. 1991). 

Data reduction was performed using the NRAO 
Astronomical Image Processing System (AIPS). At the 
first step of the data reduction, we flagged the data of 
10-channels at both edge of the band in each IF, and then 
normalized the cross-correlation by the auto-correlation. 
A standard a priori amplitude calibration was performed 
using the AIPS task APCAL based on the measurements of 
the system temperature (Tgys) by the R-Sky method during 
the observations and the aperture efficiency provided in 
the VERA Status Report in 2009. We did not employ the 
opacity correction in the task APCAL since the Tgys derived 
from the R-sky method includes the effect of atmospheric 
absorption. This amplitude calibration provides an accuracy 
of ~ 10 %, according to a number of experience of VERA 
observations (e.g., Petrov et al. 2012). Fringe fitting was 
done using the AIPS task FRJNG. After careful flagging of 
bad data and averaging 10 sec in time domain and over a 
passband, we constructed initial source model by the CLEAN 
algorithm implemented in Difmap software package and did 
self-calibration by adjusting the visibility phase to minimize 
the residual error with conserving the closure phase. After 
performing a number of iterations with the CLEAN and phase 
self-calibration and the residual error between the model and 
observed visibility converged, we switched to both phase and 
amplitude self-calibration if it works. 

For some observations, we combined multiple wv-data sets 
taken within a few days, to improve the liu-coverage. This 
is done by the AIPS task DBCON after a priori amplitude and 
phase calibrations. After combined the data sets, we performed 
the imaging as described above. 
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Source 


Alias 


z 


pc/mas * 




TOl 36+4751 


DA 55 


0.859 


7.759 


FSRQ 


J03 19+4 130 


3C 84 


0.0176 


0.344 


NLRG/Sy2 


J1230+1223 


M87 


0.004360 


0.112 


NLRG 


J 15 12-0905 


PKS 1510-089 


0.36 


5.043 


Syl/HPQ 


J1613+3412 


DA 406 


1.39712 


8.568 


FSRQ 


J1733-1304 


NRAO 530 


0.902 


7.879 


FSRQ/LPQ 


J2202+4216 


BL Lac 


0.0686 


1.289 


BLLAC 


J2253+1608 


3C 454.3 


0.859 


7.757 


FSRQ/HPQ 



* Linear size of 1 milli-arcsec. 

t FSRQ: flat spectrum radio quasars, NLRG: nan'ow line radio galaxies, Syl: Seyfert type-1, 
Sy2: Seyfert type-2, HPQ: highly polarized quasars, LPQ: low polarized quasars, BLLAC: BL 
Lac objects 



3.1. Modeling the Data 

To quantify the position, size, and flux density of core and 
jet components, we modeled the source by multiple Gaussian 
components. We fitted the visibility data with elliptical 
Gaussians that gave a good fit to the data as judged by a relative 
statistic using the modelfit task in Difmap. We sometimes 
used elliptical Gaussians or point sources if they can produce 
better fits and consistency across the epochs. Figure 2 shows 
examples of model-fit image. Note that these images are essen- 
tially similar to total intensity images produced by the CLEAN 
(Figure 3-10), demonstrating the validity of modeling proce- 
dure. 

4. Flux Density and Positional Accuracies 

It is important to check the flux density and positional accu- 
racy for the study of time variation. Among eight sources, we 
obtained a subset of two or three images observed within a very 
short timescale 1 weeks) for 3C 84, NRAO 530, 3C 454.3, 
and BL Lac. Using these four sources, we checked the consis- 
tency of flux density of each component. In Table 2 we cal- 
culated AS" = a/S for each source, where cr is a standard de- 
viation of flux density among epochs and S is the flux density 
averaged over the epochs shown in Table 3-10. The resultant 
AS* for the core component was less than 7.5%, which did not 
exceed the typical flux calibration error of VERA ^ 10% (see 
§3). On the other hand, AS of the other jet components were 
typically 20% probably because of the lack of sampling on the 
Mv-plane. Thus, we adopt 10% for the flux error of the core 
component and 20% for that of other jet components. One 
exception was that all three components (CI, C2, and C3) of 
3C 84 show very good agreements in flux density among three 
epochs (AS" < 6%). We adopt the flux error of 10% for afl 
components of 3C 84. 

We also checked the consistency of the component position 
with reference to the core, and found that the positional dif- 
ference between epochs was typically ^0.1 mas. However, 
the result of position for the component C2 of PKS 15 10-089 
was quite variable. This is probably owing to the existence of 
multiple sub-components in C2 (see §7.4). In this paper, we 
assume that typical positional accuracy is about 0.1 mas, but 
more details will be elaborated in the forthcoming paper 



Table 2. Flux error estimation 



Source 


Component 


AS-* 


3C 84^ 


CI 


0.011 




C2 


0.062 


@ 


C3 


0.057 


PKS1510-089* 


CI 


0.018 




C2 


0.215 


NRAO 530§ 


C1+C2 


0.012 




C3 


0.208 


BL Lac# 


CI 


0.075 




C2 


0.18 


3C 454.3** 


C1+C2 


0.031 



* AS = cr/S, where cr is a standard deviation of flux 
density among epochs and S is the flux density averaged 
over the epochs, t jg calculated from the data on 

2010/NOV/28, 2010/Nov/29, and 2010/Dec/04. * AS is cal- 
culated from the data on 2010/Nov/29 and 2010/Dec/04. § 
AS is calculated from the data on 2010/Nov/ll, 2010/Nov/14 
and 2010/Nov/16. # AS is calculated from the data on 
201b/Feb/8, 2011/Feb/12 and 2011/Feb/13. ** AS is cal- 
culated from die data on 2010/Nov/lO and 2010/Nov/12. 

5. Results 

5.1. Images 

In Figures 3-10, we show three images in adjacent observ- 
ing periods for each source. We convolved with an identical 
beam which was the largest beam among three epochs for each 
source. For each image, we provide the beam size, contour 
levels, and peak intensity in the caption. Image noises are gen- 
erally between 10 and 30 mJy beam^^. The noise of 3C 454.3 
image was 76 mJy beam^^ since the image dynamic range was 
limited by very bright (more than 20 Jy) core. Physical parame- 
ters of Gaussian models for these three epochs are summarized 
in Tables 3-10. All sources exhibited a bright core and a few 
additional components in the jet. Because of the lack of short 
baselines, we only detected the structure within ~ 3 mas from 
the core, but missed extended structures. 

5.2. Light curves 

In Figure 11, we plot the light curve of our GENJI pro- 
gramme for the first 7-months and the MOJAVE Ught curve 
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for the same periods. For GENJI data, we plot the light curve 
of each component described in the section 3.1 in addition to 
that of total flux density. For MOJAVE data, we plot the light 
curve of total flux density. The flux error due to the model-fit 
error is negligible as compared to the flux calibration error. We 
assumed that the flux error of MOJAVE is 5%. It is obvious 
from Figure 1 1 that the sampling interval of GENJI is much 
denser than that of MOJAVE. 

6. Consistency Check with MOJAVE on Flux Density 

In this section, we compare our GENJI light curve with the 
MOJAVE light curve. MOJAVE provides the flux density ev- 
ery few months and the observing frequency 15 GHz is nearly 
same with GENJI observing frequency 22 GHz. Therefore, the 
MOJAVE data can be a check for the fideUty of GENJI data. 

In particular for PKS 15 10-089, there are two MOJAVE data 
which were observed in nearly same periods with GENJI. The 
flux densities from GENJI are 1.92 Jy on 2010 November 29 
and 2.31 Jy on 201 1 February 26, while those from MOJAVE 
are 1.69 Jy on 2010 November 29 and 2.16 Jy on 2011 
February 27. The flux densities of GENJI are 10-20% higher 
than those of MOJAVE in both epochs. Since the most of flux 
density originates in the core and the core has a slightly in- 
verted spectrum with a spectral index (a) of 0.2 (Sokolovsky 
et al. 2010), these differences in flux density between GENJI 
and MOJAVE are reasonable. 

Another good example is BL Lac. Although we do not have 
good overlaps in observing periods between MOJAVE and 
GENJI, the trend of light curve for the 7-months is very similar 
between them. The GENJI flux density is always smaller than 
the MOAJVE flux density by ~ 1 Jy, but this difference can be 
explained by the contribution from the extended jet. A continu- 
ously jet extending up to more than 5 mas from the core is seen 
on the MOJAVE images (Lister et al. 2009). This extending jet 
cannot be time variable as compared to the core. Therefore, the 
similar trend between GENJI and MOJAVE demonstrates the 
repeatability of the flux measurement by GENJI. 

Although there are not so many overlaps, similar trend in 
Ught curve is apparent in NRAO 530 and 3C 454.3. 

In summary, our GENJI programme is promising to provide 
dense monitoring data for flux measurements. 

7. Notes on Individual Sources 

In this section, we describe the morphology and light curve 
for each source after a brief review of background information 
and previous observations. 

7.1. J0136+4751 (DA 55) 

DA 55 is classified as Flat Spectrum Radio Quasar (FSRQ) 
at z=0.859 (Healey et al. 2008). This source is listed in the 
Fenni/LAT bright 7-ray source (OFGL J0137.1H-4751) though 
it was not detected by EGRET. The 7-ray flux density of DA 55 
in 2009 mid-March reached about three times the averaged flux 
density in the first year light curve (Abdo et al. 2010). The pc- 
scale jet feature is one-sided in the VLBA 15 GHz total inten- 
sity image (Lister et al. 2009). The jet moves toward northwest 
direction. 



We detected the jet feature in the position angle ~ —45° 
within ^ 2 mas from the core (CI) (Figure 2(a) and Figure 3). 
The jet can be fitted by a single Gaussian component C2. Until 
2010, a bright jet component (a few hundreds of mJy) located 
at mas in the position angle ~ —40° from the core can be 
seen in the MOJAVE images. This outer component is not vis- 
ible in both the GENJI and the MOJAVE images because of its 
faintness since GENJI programme has been started. 

We did not obtain many data for this source so far. No sig- 
nificant change in flux density was detected. 

7.2. J0319+4130(3C84) 

The bright radio source 3C 84 (z=0.0I76: Petrosian et al. 
2007) shows a two-sided radio jets/lobes on the north and south 
side of the core in the central 5 pc (Walker et al. 2000; Asada 
et al. 2006). In mid-1980s, the radio flux density became ex- 
ceptionally bright, more than 60 Jy at centimeter wavelengths 
(e.g., O'Dea et al. 1984), and then subsequently decreased such 
that, by the early 2000s, the radio flux density decreased to 
^ 10 Jy. Since 2005, the radio flux density has started to in- 
crease in 2005 (Abdo et al. 2009a). This source was became 
to be known as a 7-ray source after the four-months operation 
of Fermi (Abdo et al. 2009a), and 3C 84 is the typical exam- 
ple of the sources exhibiting clear time variation in the GeV 
band among the misaligned 7-ray AGNs (Kataoka et al. 2010; 
Brown & Adams 201 1). MAGIC Cherenkov telescope also de- 
teced very-high-energy 7-ray emission with a spectral cut-off 
around 100 GeV (Aleksic et al. 2012). 

The central 1-pc structure consists of three components CI, 
C2, and C3 (Figure 2(b)). The images at same frequency with 
similar resolution in earlier epochs were presented by Nagai et 
al. (2010). In Figure 4, the emission slightly extending to the 
south of C3 can be seen. This extending emission became more 
visible in later epochs. An additional minor component can be 
fitted at ~ 2-mas south from the C3 in the last few epochs after 
2011 April (Figure 12). 

While there ai-e only two data by MOJAVE, GENJI provides 
the light curve in more detail (Figure 1 1(b)). There is a differ- 
ence in flux density between MOJAVE and GENJI by a few 
Jy, but this difference can be explained by the contribution 
from the extended jet that is not recovered by VERA obser- 
vations since more extended emission is detected in MOJAVE 
images (Lister et al. 2009). While the flux density increased 
steadily at 22 GHz and 43 GHz before 2009 (Nagai et al. 2010; 
Suzuki et al. 2012), the light curve presented here does not 
show flux density change within a level of 10%. The light curve 
at 37 GHz in the same periods were presented in (Nagai et al. 
2012), and the 37-GHz flux density does not change very much 
as well. 

7.3. J1230+1223(M87) 

M 87 is a giant radio galaxy at the center of the Virgo cluster 
with a well-known one-sided jet. Its proximity (z=0. 00436: 
Rines & Geller 2008) provides a linear resolution of 1 mas 
= 0.08 pc, corresponding to ~^140i?s for a mass of the cen- 
tral black hole AIbh 6 x lO^M© (Gebhardt & Thomas 
2009). This source has been routinely detected at energies of 
MeV/GeV (Abdo et al. 2009c) and TeV (e.g. Aharonian et al. 
2006; Acciari et al. 2009; Abramowski et al. 2012). During the 
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flaring activity at TeV in 2008, simultaneous VLBI observa- 
tions detected a strong flux density increase from the radio core 
(Acciari et al. 2009), whose location is specified as -^20 Rg 
from the black hole (Hada et al. 2012). On the other hand, a 
distinct multi-wavelength correlation was found for the recent 
TeV flaring event in April 2010, where the core was relatively 
stable in radio bands (Hada et al. 2012). Note that the M 87 
jet has another remarkable feature called HST-1 located at a 
deprojected distance of at least 120 pc, although its direct con- 
nection with 7-ray activity is still under hot debate (e.g. Chang 
et al. 2010; Chang et al. 2010; Giroletti et al. 2012). 

The images in Figure 5 shows a core and jet feature. 
Although we fitted the image using only one Gaussian compo- 
nent in Figure 2, it is sometimes possible to put an additional 
component in the jet. However, it highly depends on the uv- 
coverage whether we can put this component. Thus, we present 
the model-fit results only for the core (CI). 

The flux density of CI is around 1 Jy and slightly time vari- 
able with a level of 10-20% within the observing periods pre- 
sented in this paper (Figure 5). 

7.4. J] 512-0905 (PKS 1510-089) 

PKS 1510-089 is a FSRQ at z=0.36 (Thompson et al. 1990). 
This source was discovered as a 7-ray source by EGRET and 
one of the three FSRQs detected by MAGIC and H.E.S.S. 
(Cortina 2012, Wagner 2010). Since 2008 the 7-ray emis- 
sion of PKS 1510-089 became highly variable, and many rapid 
and intense flaring episodes have been detected by AGILE and 
Fermi (e.g. D' Ammando et al. 2009; D' Ammando et al. 201 1). 
In 201 1 July and October, two major 7-ray flares were detected. 
Especially during the latter flare, the source reached its highest 
7-ray flux density, becoming the second brightest AGN ever 
observed by Fermi. Orienti et al. (2012) investigate the con- 
nection of the historical 7-ray activity with lower energy bands 
and GENJI plays an important role on the paper by providing 
VLBI core light curve as a part of our first result. Radio flux 
density shows no apparent correlation with the former 7-ray 
flare, however, if the latter flare is related to radio outburst, 
the correlation indicates that the 7-ray emission region locates 
along the jet about 10 pc downstream from the central engine 
(Orienti et al. 2012). 

We detected the radio emission within ^ 3 mas from the core 
(Figures 2(d) and 6). The radio emission mostly originates in 
the core (CI), and faint jet is extending in the position angle 
^ —30°. The jet can be represented by a single Gaussian com- 
ponent C2 in our GENJI image. In the later few epochs, an 
additional component C3 can be fitted in the region between 
CI and C2 (Figure 13). C2 and C3 correspond to Nl and N2 
of higher resolution image by VLBA (Fig.4 of Orienti et al. 
201 1), respectively. 

We do not include the flux density of C3 in Figure 11(d) 
since C3 appeared to be visible only in later few epochs. Total 
flux density had decreased from 2.22 Jy on 2010 November 
1 to 1.58 Jy on 2010 December 14 (Figure 11(d)). It became 
to be bright in 2011 February (2.04 Jy on 2011 February 11 
and 2.26 Jy on 2011 February 26), and then subsequently de- 
creased to about 1 .5 Jy in 201 1 April. Overall trend is similar to 
the MOJAVE, but our GENJI light curve provides more detail 
change in flux density. 



7.5. J1613+3412 (DA 406) 

DA 406, which is classified as FSRQ, is the highest-redshift 
object in the GENJI source list (z = 1.397: Schneider et al. 
2007). One-sided jet is extending to the southern-east direc- 
tion in the VLBA 15 GHz image (Lister et al. 2009). A bright 
component is located at about 4.5 mas from the core. This 
source was detected with both EGRET (3EG J1614H-3424) and 
Fenni/LAT (IFGL J1613.5H-341 1). The 7-ray flux density of 
this source in EGRET-era was clearly brighter than in Fermi- 
era (Finer & Kingham 1997; Abdo et al. 2010). 

Our GENJI images show a bright radio core and one-sided 
jet extending in the position angle of ^ 170° (Figures 2(e) and 
7). The central ~ 2 mas structure can be fitted with two com- 
ponents CI and C2. Another component C3 can be fitted in the 
jet feature at 4.5 mas from the core. No significant motion of 
these components can be seen by GENJI 7-month monitor. 

Our GENJI light curve shows that the flux densities in the 
later periods (from 2011 April 22 to 2011 May 25) clearly 
increased from those in earlier epochs (Figure 11(e)). The 
MOJAVE provides only one data point in these periods, and 
therefore this change in flux density is not detected. 

7.6. J1733-1304(NRAO 530) 

NRAO 530 (also known as PKS 1730-13) is a wefl-known 
blazar at z=0.902 (Healey et al. 2008) classified as a FSRQ 
source. NRAO 530 has a strong variability at various wave- 
length such as radio (Bower et al. 1997), optical (Webb et 
al. 1988), X-ray (Foschini et al. 2006) and the 7-ray bands 
(Mukherjee et al. 1997; D' Ammando & Vandenbroucke 2010). 
In 7-ray band, NRAO 530 was identified with the EGRET 
source 2EG 1735-1312 (Thompson et al. 1995). After the 
launch of Fermi, it did not appear as a bright source (Abdo et 
al. 2009b) but as a relatively quiescent source IFGL J 1733.0- 
1308 (Abdo et al. 2010b). However, its 7-ray activity in- 
creased at the latter half of 2010 and NRAO 530 flared be- 
tween October 31 and November 2, 2010 (D' Ammando & 
Vandenbroucke 2010). On pc-scale, NRAO 530 has a core-jet 
structure extending to the north from the core with a wiggled 
trajectory on a scale of about 30 mas (Lu et al. 201 1). 

We detected a structure within 5 mas from the core (Figures 
2(f) and 8). The observed structure is consistent with previous 
VLBI observations. It mainly consists of two major compo- 
nents CI, C2 and two faint jet components C3 and C4 (Figure 
2). 

In Figure 1 1, we show the sum of CI and C2, since the dis- 
tance between CI and C2 is very small and the flux density 
of CI can be correlated with that of C2. In some epochs, we 
could not distinguish C2 from CI or we could not detect C3 
and C4 presumably owing to a lack of uu-coverage. We do 
not include the light curve of C4 in Figure 1 1 since the fit is 
not quite robust. Our observations started 9 days after the 2010 
October 7-ray flare. Our results show that the core flux density 
of NRAO 530 had increased from 2.7 Jy to 4.0 Jy in the 
end of 2010. The flared core flux density had been relatively 
stable until 201 1 May and started to decrease in 201 1 May. The 
change in flux density in this period is difficult to follow only 
with MOJAVE. NRAO 530 is one of the most successful exam- 
ples to show how intensive monitoring provided by GENJI is 
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important to study the connection between 7-ray flare and the 
change in radio flux density. 

7. 7. J2202+4216 (BL Lac) 

BL Lacertae (z=0.0686: Sazonov et al. 2007) is a prototypi- 
cal blazar which has been well studied across the entire electro- 
magnetic spectrum. This source has been detected by EGRET 
several times and now routinely detected by Fermi (Abdo et 
al. 2011). Above 200 GeV, a significant detection was made 
by MAGIC (Albert et al. 2007). VLBI structure of this source 
is characterized by the bright core and the extended jet propa- 
gating southward (e.g. Jorstad et al. 2005; Bach et al. 2006). 
Recent multi-wavelength study detected a multiple 7-ray flar- 
ing event which is followed by a emergence of a superluminal 
knot from the core, together with a continuous rotation of the 
optical polarization angle (Marscher et al. 2008). These lead to 
the scenario that the flaring event in this source is produced at 
large distance (more than 1 pc or 10^ Rs) downstream of the 
central engine, although such an interpretation is still contro- 
versial yet (Tavecchio et al. 2010). 

We detected the bright core (CI) and a jet feature extend- 
ing to a few milli-arcsec from the core. The jet feature can be 
represented as a single Gaussian component C2 (Figure 2(g)). 

The difference in flux density between our GENJI light curve 
and MOJAVE light curve (see Figure 11(g)) can be explained 
by the contribution from the extended jet since a continuous 
jet extending up to more than 5 mas from the core is seen on 
the MOJAVE images (Lister et al. 2009). The flux density de- 
creased from 4.0 Jy on 2010 November 22 to 2.5 Jy on 2011 
January 15, and then small rise and decay in flux density is seen 
around 201 1 February to March. In later epochs, the flux den- 
sity increased slightly from 2.6 Jy on 2011 March 2 to 3.2 Jy 
on 2011 May 22. 

7.8. 3C 454.3 

A weU-known FSRQ source 3C 454.3 at z= 0.859 (Sazonov 
et al. 2007) is one of the brightest extragalactic radio sources. 
After 2000, the radio flux density of 3C 454.3 had been quite 
time variable and showed increased activity many times. In 
particular since 2005, a remarkable flaring activity was seen in 
this source across the whole electromagnetic spectrum from ra- 
dio to 7-rays. During the last four years, this source has exhib- 
ited more than one 7-ray flare every year, becoming the most 
active 7-ray blazar in the sky. The largest flare occurred in 
November 2010. In this flare, 3C 454.3 reached a peak flux 
density more than a factor of 6 higher than the flux density of 
the brightest steady 7-ray source, the Vela pulsar, and more 
than a factor of 3 brighter than its previous super-flare on 2- 
3 December 2009 (Abdo et al. 2011; Verceflone et al. 2011). 
On pc-scale, 3C 454.3 has a core-jet structures extending to 
the west. Previous VLBI observations showed jet components 
propagating along different trajectories (Lister et al. 2009). 

We detected a structure within 1 mas from the core. It mainly 
consists of two major components CI, C2 (Figure 2(h)). 

In Figure 10, we show the sum of flux densities of CI and 
C2, since the distance between CI and C2 is very small and 
flux density of CI can be correlated with C3. Our observations 
started just before the largest 7-ray flare in November 2010. 
Our results show that the core flux of 3C 454.3 started to in- 



crease in 2010 December and peaked around 2011 February. 
This is consistent with the idea that the 7-ray flaring region is 
located at the upstream of the radio core and the flux density 
increase in the radio core is the result of the propagation down 
of the flaring component. If we assume a jet speed of F = 10 
(F is the bulk Lorentz factor), the 7-ray flaring zone should be 
~ 0.05 pc upstream of the 22-GHz radio core. We note that this 
rise and decay in flux density is visible from the GENJI light 
curve but is not very apparent from the MOJAVE light curve 
(Figure 11). 

8. Summary 

We presented the first 7-month results from the GENJI pro- 
gramme. In this paper, we particularly focused on the fi- 
delity of the image and flux density measurement by compar- 
ing nearby epoch data, and we also presented the consistency 
with the MOJAVE. We will discuss the detail of light curve and 
kinematics for individual sources in separate papers. In partic- 
ular, we successfully obtained the data right before/after the 
7-ray flares reported in PKS 15 10-089, NRAO 530, 3C 454.3. 
The 22-GHz data from the GENJI are unique and they provide 
complementary information with other monitoring projects in 
radio bands. 
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Fig. 1. Examples of uv-covQiage (a) DA 55, (b) 3C 84, (c) M 87, (d) PKS 1510-089, (e) DA 406, (f) NRAO 530, (g) BL Lac, (h) 3C 454.3. 
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Fig. 2. Examples of Gaussian model fit of (a) DA 55, (b) 3C 84, (c) M 87. (d) PKS 1510-091, (e) DA 406, (f) NRAO 530, (g) BL Lac, and (h) 3C 454.3. 
Total intensity image is created by CLEAN and model-fit image is produced by restoring the model-fit components. 



Nagai et al. 



9 



DA55 at 22.206 GHz 2011 Apr 22 




DA55 at 22.202 GHz 2011 May 06 




DA55 at 22.202 GHz 2011 May U 




-2 
Right Ascension (mos) 



Right Ascension (m 



-2 
Right Ascension (mosj 



Fig. 3. Total intensity images of DA 55. All images are convolved with an identical synthesized beam of (1.2x0.8) mas at the position angle of —50 deg. 
The contours are plotted at the level of 0.025 x{-l, 1, 2, 4,. . ., 64) Jy beam~^. The peak intensities of three images are 1.80, 1.83, and 1.84 Jy beam~^, 
respectively. 
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Fig. 4. Total intensity images of 3C 84. All images are convolved with an identical synthesized beam of (1.26 X 0.861) mas at the position angle of 
—67.4°. The contours are plotted at the level of 18.1 X (—1, 1, 2, 4, ... , 32) mjy beam^^. The peak intensities of three images are 7.05, 6.79, and 
7.27 Jy beam~^, respectively. 




Fig. 5. Total intensity images of M 87. All images are convolved with an identical synthesized beam of (1.36 X 0.87) mas at the position angle of —43°. 
The contours are plotted at the level of 18 X (— 1, 1, 2, 4, . . . , 32) mJy beam~^. The peak intensities of three images are 0.95, 0.95, and 0.95 Jy beam~^, 
respectively. 
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Fig. 9. Total intensity images of BL Lac. All images are convolved with an identical synthesized beam of (1.3 X 1.0) mas at the position angle of —40°. 
The contours are plotted at the level of 30 X (—1, 1, 2, 4, ... , 64) mJy beam^^. The peak intensities of three images are 2.30, 2.08, and 2.30 Jy beam~^, 
respectively. 




Table 3. Summary of Gaussian model parameters of DA 55. 



Epoch 


Component 


Flux [Jy] 


Radius [mas] * 


B [deg] T 


Major [mas] ^ 


Axial ratio^ 


[deg] ** 


2011/Apr/22 


CI 


1.803 


0.012 


115.4 


0.212 


0.000 


-79.4 




C2 


0.161 


0.971 


-49.8 


0.606 


1.000 


-21.8 


2011/May/06+May/09 


CI 


1.820 


0.012 


123.6 


0.107 


0.000 


71.4 




C2 


0.162 


0.992 


-52.3 


0.867 


1.000 


17.1 


2011/May/19 


CI 


1.877 


0.012 


123.0 


0.291 


0.140 


-72.0 




C2 


0.197 


1.050 


-48.8 


0.726 


1.000 


0.0 



* The radial distance of the component center from the center of the map. t The position angle of the center of the component, -f The FWHM major 
axis of the component. § The ratio of the minor axis to the major axis. ** The position angle of the major axis. 
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Table 4. Summary of Gaussian model parameters of 3C 84. 



[Vol. , 



Epoch 


Component 


Flux [Jy] 


Radius [mas] 


[deg] 


Major [mas] 


Axial ratio 


4> [deg] 


2010/NOV/28 


CI 


5.493 


1.409 


-4.6 


0.744 


0.451 


-20.7 




C2 


4.956 


0.595 


-75.8 


0.761 


1.000 


-101.3 




C3 


6.895 


0.165 


108.9 


0.680 


0.308 


12.7 


2OIO/N0V/29 


CI 


5.447 


1.449 


-5.8 


0.793 


0.423 


-19.6 




C2 


5.020 


0.602 


-73.1 


0.757 


1.000 


-90.0 




C3 


6.166 


0.140 


99.3 


0.652 


0.000 


14.7 


2010/Dec/4 


CI 


5.568 


1.361 


-3.4 


0.768 


0.401 


-21.3 




C2 


5.537 


0.568 


-77.4 


0.867 


1.000 


-104.0 




C3 


6.690 


0.191 


119.5 


0.595 


0.382 


19.3 


Table 5. Summary of Gaussian model parameters of M 87. 


Epoch 


Component 


Flux [Jy] 


Radius [mas] 


[deg] 


Major [mas] 


Axial ratio 


[deg] 


2OIO/N0V/O4 


CI 


1.131 


0.046 


-77.3 


0.796 


0.384 


-60.8 


2010/Dec/04 


CI 


1.120 


0.035 


-85.0 


0.695 


0.515 


-59.5 


2010/Dec/16 


CI 


1.205 


0.031 


-73.9 


0.716 


0.694 


-52.0 


Table 6. Summary of Gaussian model parameters of PKS 1510-089. 


Epoch 


Component 


Flux [Jy] 


Radius [mas] 


[deg] 


Major [mas] 


Axial ratio 


[deg] 


2OIO/N0V/I8 


CI 


1.968 


0.008 


-88.2 


0.538 


0.349 


-33.6 




C2 


0.129 


2.523 


-30.8 


1.614 


0.292 


-17.3 


2OIO/N0V/29 


CI 


1.817 


0.014 


-47.1 


0.618 


0.368 


-30.7 




C2 


0.112 


2.474 


-36.1 


0.777 


0.546 


-25.6 


2010/Dec/04 


CI 


1.771 


0.022 


-50.7 


0.544 


0.436 


-30.4 




C2 


0.081 


2.118 


-33.6 


0.000 


1.000 


0.0 



Table 7. Summary of Gaussian model parameters of DA 406. 



Epoch Component Flux [Jy] Radius [mas] 6 [deg] Major [mas] Axial ratio [deg] 



2010/Dec/21 


CI 


1.256 


0.007 


7.647 


1.491E-08 


1.000 


-3.5 




C2 


0.110 


4.355 


157.1 


0.826 


1.000 


8.5 




C3 


0.041 


1.043 


170.5 


0.607 


1.000 


-12.2 


2011/Jan/09H-Jan/15 


CI 


1.113 


0.034 


-2.8 


1.761e-07 


1.000 


-8.5 




C2 


0.049 


4.303 


158.8 


0.194 


1.000 


0.0 




C3 


0.149 


0.526 


177.7 


0.533 


1.000 


-8.5 


2011/Feb/02 


CI 


1.281 


0.012 


-2.1 


0.130 


1.000 


-36.2 




C2 


0.056 


4.202 


160.1 


0.660 


1.000 


-10.1 




C3 


0.071 


0.798 


168.8 


5.067e-07 


1.000 


-5.6 
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Fig. 11. The first 6-month Ught curve of (a) DA 55, (b) 3C 84, (c) M 87, (d) PKS 1510-089, (e) DA 406, (t) NRAO 530, (g) BL Lac, and (h) 3C 454.3. 
For the comparison, MOJAVE flux data is also plotted. 
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Table 8. Summary of Gaussian model parameters of NRAO 530. 



[Vol. , 



Epoch 


Component 


Flux [Jy] 


Radius [mas] 


9 [deg] 


Major [mas] 


Axial ratio 


J. r J 1 

9 [deg] 


2010/Nov/ll 


CI 


2.733 


0.118 


-170.1 


0.000 


1.000 


0.0 






i .4Uj 




0.0 


U.UUU 


i .UUU 


U.U 




C3 


0.118 


2.682 


-3.9 


1.723 


0.000 


-36.2 


2010/NOV/14 


CI 


2.594 


0.129 


-177.1 


0.000 


1.000 


0.0 




C2 


1.573 


0.228 


2.7 


0.000 


1.000 


0.0 




C3 


0.174 


2.715 


-3.8 


1.517 


0.516 


-30.886 


2010/NOV/16 


CI 


2.925 


0.092 


-168.0 


0.000 


1.000 


0.0 




C2 


1.142 


0.262 


11.9 


0.000 


1.000 


0.0 




C3 


0.130 


2.639 


-6.1 


1.869 


0.456 


-46.8 



Table 9. Summary of Gaussian model parameters of BL Lac. 



Epoch 


Component 


Flux [Jy] 


Radius [mas] 


9 [deg] 


Major [mas] 


Axial ratio 


(f) [deg] 


2011/Feb/08 


CI 


2.477 


0.054 


171.7 


0.598 


0.000 


41.5 




C2 


0.593 


2.125 


-165.5 


0.788 


0.743 


-6.9 


2011/Feb/12 


CI 


2.175 


0.006 


171.6 


0.330 


0.000 


45.4 




C2 


0.843 


2.038 


-167.8 


2.137 


0.306 


10.3 


2011/Feb/13 


CI 


2.489 


0.015 


158.1 


0.423 


0.333 


20.0 




C2 


0.680 


1.945 


-169.0 


1.265 


0.445 


17.2 
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Table 10. Summary of Gaussian model parameters of 3C 454.3. 



Epoch 


Component 


Flux [Jy] 


Radius [mas] 


6 [deg] 


Major [mas] Axial ratio 


4> [deg] 


2010/Nov/lO 


CI 


21.039 


0.053 


79.1 


1 







C2 


3.987 


0.363 


-102.7 


1 





2OIO/N0V/I2 


CI 


21.230 


0.041 


80.9 


1 







C2 


2.771 


0.423 


-100.2 


1 





2010/Dec/ll 


CI 


21.487 


0.046 


76.0 


1 







C2 


3.456 


0.390 


-104.3 


1 







Fig. 12. Modelfit image of 3C 84 on 20 11 April 21. A point source 
can be fitted at ~ 2-mas south from the C3. 




Fig. 13. Modelfit image of PKS 15 10-089 on 2011 May 12. 



